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SUMMARY 


Previous simulator studies have shown that a proposed supersonic transport 

sideslin e * hihit ? undes l ra ble lateral motions during landing approach. Large adverse 
desiip excursions and large peak lateral acceleration at the pilot’s station 
occurred during rolling maneuvers of the unaugmented airplane. station 

that nrnv 1 !! 8 ®, tUd /* C °“ tro1 theory has been applied to determine feedback gains 

to aillron ^ defle^tTon i 8t * blll H ty characteristics and satisfactory transient response 
^ dafleatlon in P ut ‘ However, the peak value of lateral acceleration at the 
pilot station does not satisfy a proposed criterion during a rolling maneuver. 

1 fl t P r!i Ptlnial i regulator theory does not provide a significant reduction of the peak 
extrpm! acceleration • The weighting matrices in the performance index were given an 
variation, but the effect on lateral acceleration remained insignificant 
Jj 0 o /se °° P control provided the desired bank angle (30°) with the desired roll rate 

larg^ adlersf sKeI?i a satisfactory level of lateral acceleration. However, a 
ntliti dve / se sldesl± p angle was required. In addition, the yawing velocity was 

maneuver . SeC ° ndS lat6ral acce l e ration changed sign during the 

! be P rob l em Persists and, as suggested in previous studies, perhaps relaxed 
age rolling S velocTtv P ° Sed 1* ^ C ^ lteria cou l d allow any combination of a lower aver- 
lateral Selection ^ at ^ “« U »•* 


INTRODUCTION 


jssS 

_ * h \ f 1 . tion characteristics were considerably improved. In response to a 

evident' A £dV« TJwJf 186 ^ lateral accelerations at the pilot station were 
the peak lateral** 1 d stability and control augmentation system successfully reduced 

2en a .tip Teel* Z„t ' ^ r ° UinS Velocity Was also considerably reduced 

^ ® el input w as applied to the airplane. Another simulation (Ref 2^ of 

to extreme i rP f an6 i With 1 Stabillty au 8 menta tlon systems attributed high pilot ratings 
inv mane & eral acce l er ations that were produced at the pilot station during roll- 
ing maneuvers. g Suggested ways to alleviate the undesirable "peak lateral acceSratton 

restrict rolling maneuvers; a'nd all “ some* adfersfsideTlip!^ lat6ral aCCeleratl ° n; 

control 8 toeoVv 7 system has bee " designed analytically by use of modal 

y. The requirements of reference 3 are satisfied with use of the 
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system. Optimal regulator theory does not provide significant reduction of lateral 
acceleration at the pilot's station. Open-loop control was applied, therefore, to 
find control inputs that could satisfy the criteria (Refs. 2 and 4) for peak lateral 
acceleration. 

The airplane was assumed to be rigid and control inputs were aileron and rudder 
deflection angles . 


SYMBOLS 


Values are given in the International System of Units (SI). Dots over symbols 
denote differentiation with respect to time. All calculations are based on the 
airplane body axes. 


^ a y^PS 

(a ) 
y eg 

B 

C 


n 

Si 

c o 

U 6a 

D 

F,S,G 

g 

j 

K 

M 

I 

n-m 

I X> I Y * I z 

x xz 

P>r 

r 2 ,r 4 

t 

T 2 

T l/2 


V 

V 


matrix of coefficients in equation (1) 

matrix of augmented coefficients as in equation (A-9) 
lateral acceleration at the pilot's station, g units 

lateral acceleration at the center-of-gravity of the airplane, 
g units 

control effectiveness matrix in equation (1) 
matrix of coefficients in equation (2) 

mean aerodynamic chord, m 
side force coefficient 

rolling moment coefficient 

yawing moment coefficient 

matrix used in equation (A-8), (i=l,2,3,4) 
coefficient defined in table II 

matrix of coefficients in equation (2) 
matrices defined by equation (A-7) 
free fall acceleration, m/sec^ 

/-[ 

feedback gain matrix 

matrix with elements m i;l , m ik , m^, l<i<n, l<k<n 
unity matrix having dimension (n-m)x(n-m) 

moment of inertia about X, Y, and Z body axes, respectively, kg— m^ 
product of inertia, kg-ra^ 

rolling and yawing velocities, respectively, deg/sec 

real space of 2 and 4 dimensions, respectively 
time, sec 

time to double amplitude, sec 
time to half amplitude, sec 
control vector 

control vector of commanded inputs 
modal matrix 

airplane airspeed, m/sec 
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V 11 ,V 22 ,V 23 » 
v 32 ,v 33’ v 44 
V R ,V S 

V DR,R ,V DR,I 

X 

x(t),x(o) 



characteristic vectors (i=l,2,3,4) 
elements of modal matrix in equation (8) 

characteristic vectors associated with the roll mode and spiral 
mode, respectively 

real and imaginary vectors, respectively, of the characteristic 
vector associated with the Dutch-roll mode 
state vector defined after equation (1) 

state vector as a function of time and initial value of the state 
vector, respectively 

longitudinal and vertical distances, respectively, from the airplane 
center-of-gravity to pilot station, m 
component vectors of a characteristic vector (i=l,2,3,4) 

angle of attack, deg 
angle of sideslip, deg 

aileron deflection, positive for right roll command, deg 
flaperon deflection, deg 

aileron deflection caused by pilot action, deg 
trailing-edge flap deflection, deg 
rudder deflection, deg 

rudder deflection caused by pilot action, deg 

deflection of spoiler-slot and inverted spoiler-slot deflectors, deg 

elements of modal matrix in equation (8) 

Dutch roll mode damping ratio 

characteristic values (i=l,2,3,4) 

characteristic value associated with roll mode 

characteristic value associated with spiral mode 

matrix of characteristic values defined in equation (A-10) 
elements of modal matrix in equation (8) 

roll mode time constant, sec 

undamped natural frequency of Dutch roll mode, rad/sec 


PRELIMINARIES 


Mathematical Model 

The airplane considered in this report is the same as the baseline concept of 
reference 1. A sketch of the airplane is shown in figure 1. The mass and 

dimensional characteristics and aerodynamic characteristics (a«8°) are given in 
tables I and II, respectively. The linearized equations of lateral motion for a 

f pp 5°? ch con ^8uration («fi aps =40 o , V 0 =78.71 m/sec and Altitude-91 .44 

meters), in body-axes, are given by the following equation. 
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x = Ax + Bu 


(i) 


where 



p 


<5 


r 


a 

X = 

6 

; u = 

6 


♦ 


r 


-.5676 

.8599 

-2.193 

0 


.8295 

.2106 

-.0069 

-.1475 

.2301 

0 

; B = 

-.0135 

-.1159 

.1950 

-.9706 

-.0799 

.1234 

-.0065 

.0111 

1.0 

.1405 

0 

0 


0 

0 


Additionally, the lateral acceleration at the pilot's station is given, in g units, 
by 

( a y ^PS = CX + Du (2) 

where 


C - [.1405 -.0879 -.6732 0] 


D = [.2915 -.3311] 

The pilot's station is located 44.2 meters forward (x-direction) of the center-of- 
gravity (eg) and 4.78 meters above (z-direction) the eg. The effect of pilot dis- 
placement in the y-direction and the effect of pitch rate on lateral acceleration at 
the pilot's station have been neglected. In terms of distance from the eg to the 
pilot's station, the lateral acceleration is, with the above assumptions. 


• a y ^PS 


= (a ) + 

y eg 


rx - pz 


(3) 


where x » 44.2 m, and z = -4.78 m. 


Lateral-Directional Characteristics (Unaugmented Airplane) 

The characteristic values for the supersonic transport airplane modelled in the 
preceding section are the following (See also table III): 

Roll Mode: 1 ■ -.611 

R 


Dutch Roll Mode: A ■ -.077 ± .821 j 

DR 

Spiral Mode: A - -.031 

u 
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The roll mode time constant, t r =-1/X r « 1 .64 sec, exceeds the requirement of refer- 
ence 3 (See Table III). Damping of the Dutch-roll mode, moreover, is too light for 
this airplane. The effect of these characteristics is shown in the response of the 
airplane to a moderate aileron angle input (6 a =15°) (Fig. 2). Since the rolling 
velocity, p, is oscillatory, the Dutch-roll mode and roll mode are coupled, and the 
rolling velocity does not approach a steady-state value in a desired manner. The 
spiral mode, although stable, shows an undesirable lag in the response of the yawing 
velocity. Table III compares the airplane characteristics with those specified in 
reference 3. 

At the pilot’s station, the calculated lateral acceleration is quite large for 
the moderate control input (See Fig. 2(b)). The criterion given in references 2 and 
4 would certainly be unsatisfied during the prescribed rolling maneuver. The crite- 
rion of reference 4 states in part "Lateral acceleration at the pilot station shall 
not exceed a level of ±0.075 g peak, and the critical passenger station shall not 
exceed ±0.05 g peak. These levels shall be met for all normal maneuvers including 
30° bank and capture using an average roll rate of 5°/sec in cruise and 10°/sec at 
landing. If unpiloted time studies are conducted, the wheel input should be a 0.5 
second ramp of magnitude sufficient to produce the specified average roll rates." In 
reference 2, a criterion was suggested as follows: "A lateral acceleration criterion 
for large transports in the landing approach was developed. The parameter 
Ny p iiot /^nax correlates well with pilot rating. A value less than .012 g' s/deg/ 
sec was required for satisfactory ratings and .035 g f s/deg/sec for acceptable 
ratings." The largest part of the lateral acceleration comes from the second term of 
equation (3) which depends on the location of the pilot with respect to the eg. The 
airplane characteristics that have been mentioned are essentially those given in ref- 
erence 1 for the unaugmented airplane. Differences arise primarily from slight dif- 
ferences in modelling. The simulations of reference 1 lead to high pilot ratings. 
As stated, the major objections were (1) unacceptable large adverse sideslip excur- 
sions in turns; (2) easily excited, lightly damped Dutch roll mode; (3) poor roll and 
heading control; and (4) sluggish roll response with low roll damping. Also men- 
tioned in reference 1 was that large adverse sideslip excursions occurred during 
rolling maneuvers. 


STATEMENT OF THE PROBLEM 


From the preceding discussion, it is evident that a stability augmentation 
system is needed for the airplane. Also, a controller is needed to satisfy the cri- 
terion for lateral acceleration level during a rolling maneuver. In the simulations 
of references 1 and 2, several augmentation systems were used to satisfy the flying 
quality requirements in reference 3. However, although the lateral acceleration 
level was not sufficiently reduced to produce satisfactory handling qualities, it 
should be noted that the introduction of a first-order lag to the roll-rate command 
signal did produce acceptable handling qualities (6.5 > pilot rating > 3.5). 

The objective of the present report is: (1) apply modal control theory in the 

design of SAS; and (2) apply optimal regulator theory to the combined airplane/SAS 
system and, thereby, obtain theoretical control inputs to augment stability and to 
alleviate lateral acceleration levels. This two step process was followed because 
modal control theory provides a means of calculating feedback gains that result in 
pole placement and partial mode separation. Optimal regulator theory provides the 
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feedback gains that result in lower acceleration levels. Because both methods are 
regulator design methods, ultimately the gains from each method are effectively added 
in the feedback loop. Consequently, the optimal regulator will tend to negate the 
effect of the modal control regulator. If applied in the reverse order, the modal 
control regulator will tend to negate the effect of the optimal regulator. The 
overall result of the present application was considered best. 


STABILITY AUGMENTATION 


Analytical expressions for the gains that will produce any desired character- 
istic values of an augmented system are well known (Refs. 5, 6 and 7). The feedback 
gains are obtained as solutions of linear algebraic equations if the linear dynamical 
system has the property of being completely controllable using a single linear con- 
trol function. When more than one linear control function is available, additional 
system characteristics can be obtained. Response of the dynamical system can be 
shaped by partial specification of the characteristic vectors of the augmented system 
(See the Appendix for an outline of the theory of Ref. 8). How much shaping of the 
responses can be obtained depends on the number of independent control inputs . 
Specifically, if n is the dimension of the state space and m is the dimension of the 
control space, nm elements of the matrix of characteristic vectors of the closed-loop 
system are arbitrary. As a design procedure, these elements may be specified. Con- 
sequently, linear algebraic equations are available to calculate feedback gains that 
assure that the dynamical system will have desired stability and response 
characteristics . 


Characteristic Values 

The choice of characteristic values must satisfy the criteria summarized in 
table III. Since the required stability parameters are given as inequalities, the 
choices come from a wide range of values. The spiral mode of the airplane is stable 
and, of course, satisfies the criterion of reference 3. Therefore, the character- 
istic value associated with this mode, Xg=-.031, need not be changed. 

On the other hand, the time constant of the roll mode does not satisfy the 
appropriate criterion. The largest acceptable value of the roll mode time constant 
is t r <1.4 with the corresponding value of the characteristic value of X R <-.714. 
Obviously, a value of X R «-.75 should be satisfactory. For a 15° commanded aileron 
input (See Fig. 3.), the roll performance, as expressed by p and <|>, is almost halved 
when X R is decreased from -.75 to -1.5. Also, (ay)pg is almost halved. How- 

ever, a choice of X R =-1.5 (t r =>.67) was made so that the corresponding character- 
istic vector would be more acceptable. 

The Dutch-roll mode only needs more damping (That is, ^d^nd needs to be 
doubled, at least). Consesquently, Cd“*182, u> n< j , **825, Sd^nd 0 *^ satisfies 

the criterion of reference 3. Figure 4 shows that the level of lateral acceleration 
at the pilot '8 station decreases with decreasing undamped frequency, So, 

?d“»351, <%d"*^27, and Cd^nd"*^ were chosen for the Dutch-roll mode. The 

corresponding characteristic value is Xp R ».15±0.4j. Although feedback gains 
required to effect this change in damped frequency might be large, consideration of 
the lower lateral acceleration was considered to be overriding. 
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In summary, the following characteristic values were chosen as satisfactory. 


DR 


-1.5 (t r - 

.67) 


-.031 (T 1/2 

“ 22.4 sec.) 


- -.15 ± .4j 

? d = 

.351 


0) , ** 

.427 


nd 



'd"nd ' 

.15 


Characteristic Vectors 


Some definitions, and the theory necessary for understanding the following 
discussion, are given in the Appendix (see also Refs. 8 and 9). The solution of the 
closed-loop system 


. x = (A+BK)x + Bu c 

In terms of the modal matrix V can be expressed as 

x(t) = V exp (At)V ^x(O) 

when u c =0 (See Appendix for definition of A). If x is defined as 

x = V _1 x(0), 


then 


x(t) = V exp(At)x 


(A) 

(5) 

( 6 ) 
(7) 


Although the entries of the modal matrix occur nonlinearly in the x terms, dominance 


arguments show that, if V for 

the airplane 

is the 

form 




v n 

£ 1 

e 2 

e 3 



V 


s i 

V 22 

V 23 

e 4 


(8) 



S 2 

V 32 

V 33 

£ 5 




_ 6 3 


«5 

V 44 




whenever t and <5 values are small compared with v values, then V - *- will be dominant 
and cross-coupling between modes will be minimal. Note that an nxn matrix, M, is a 
dominant matrix (Ref. 10) if 


| m ii| ^ ^ | m ik|* 1 .1 1 ^ n ( row dominance) (9) 

k*i 
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and 



k*i 


< i < n (column dominance) 


( 10 ) 


are satisfied. (| | means absolute value.) 
The modal matrix for the airplane is 


or 



.5209 

.6163 

.0 

.0469 

TT 

.0023 

-.0636 

.0453 

.1137 

V — 

.0111 

-.1629 

-.2130 

.0561 


-.8535 

-.0606 

-.7338 

.9908 


< 

II 

T* 

5* 

V DR,R 

V DR,I V S ] 



(ID 


( 12 ) 


where vr is the characteristic vector of the roll mode, vdr r and vdr i are the 

real and imaginary parts, respectively, of the characteristic vector of the Dutch- 

roll mode, and vg is the characteristic vector of the spiral mode. The character- 
isti c vect ors are all assumed to be normalized by dividing each element of the vector 

by / v^*Vi where v-£* is the complex conjugate of v^. 


Also, the corresponding matrix of characteristic values for the airplane is 




-.611 

0 

0 

0 



A 


0 

-.077 

,821 

0 


(13) 


0 

-.821 

,077 

0 




0 

0 

0 

.031 



where 


V 

II 

-.611 







V- 

-.077 ± j 

(.821) 



(14) 



x s " 

-.031 





The large entry 

In 

V DR,R 

indicates 

that 

the 

Dutch-roll mode 

interacts 


strongly with the roll mode. Figure 2 shows this Interaction as previously men- 
tioned. Not so evident, however, is the strong interaction of vdr j and vr on 
the spiral mode. 


Specification of the i=l,2,3,4, vectors for any mode must satisfy equation 

A-6 for a given characteristic value, X^. For the roll mode, Xr»-.72 is a satis- 
factory characteristic value. However, the corresponding characteristic vector is 
not satisfactory because it does not satisfy equation (10). Consequently, from cal- 
culation of w-£ for many combinations of X^ and z^, the characteristic value for 
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the roll mode was selected as X R =-1.5 with the corresponding characteristic vector 
of 


.8300 

.0 

-.0705 

-.5533 


The Dutch-roll mode characteristic vector, for Xjj R *-.15±.4j, is 


(15) 


DR,R 


.0 

.2691 

-.5800 

.0518 


and 


DR, I 


.0 

.2691 

.7094 

-.1140 


(16) 


The selected characteristic value, X s =-.031, yields a corresponding characteristic 
vector of 


V S = 


-.0458 

.1093 

.1376 

.9834 


(17) 


Of course, many other characteristic vectors are available for Xg=-.031. The vg 
given, however, resulted in a more favorable feedback gain matrix than did most other 
values of vg. The gain matrix resulted in small rudder angles initially in 
response to aileron angle input. 


is 


the 

selected 

characteristic 

values 


.8300 

.0 

.0 

-.0458 


.0 

.2691 

.2691 

.1093 

V = 

-.0705 

-.5800 

.7094 

.1376 


-.5533 

.0518 

-.1140 

.9834 


(18) 


The modal matrix, V, exhibits the required dominance and the corresponding matrix of 
characteristic values, A, satisfies the criteria of reference 3. The matrix of char- 
acteristic values may be written as 


-1.5 .0 



.0 .0 


The resulting feedback gain matrix is 


-1.176 -1.172 

-.016 .508 


.0 .0 

.40 .0 

-.15 .0 

.0 -.03 


2.569 -.097 

.252 .055 


(19) 


( 20 ) 


Response of the augmented airplane to an aileron angle input of 15° is presented in 
figure 5. Comparison with the response for the unaugmented airplane given in fig- 
ure 2, shows that a more desired response has been achieved. However, the peak 
lateral acceleration at the pilot's station is still too large. 
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LATERAL ACCELERATION ALLEVIATION 


The airplane, augmented as described in the preceding section is represented by 

x = (A+BK)x + Bu 


or 


with 


x « Ax + Bu 


A = A + BK 


( 21 ) 


( 22 ) 


The controlled variable is lateral acceleration at the pilot's station and is repre- 
sented as 


(ay) ps = Cx + Du 


(See Eq. (2) and (3).) 

A control law, given by 

minimizes a performance index 


-IT T 
u = R [B P + S ]x 


J = Iq ( a y)p S + u T RjU)dt 


or 


J = /q (x^Qx + 2x X Su + u X Ru)dt 
where the matrix P satisfies the matrix Riccati equation 

0 = -PA - A T P - Q + (PB + S)(B T P + S T ) 


In equation (24) and equation (25), 


T 

Q -■ C C 


T 

S = CD 


(23) 


(24) 


(25) 


(26) 


(27) 


R « R x + D D 
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Lateral acceleration at the pilot’s station of the supersonic transport airplane 
is much too large during a rolling maneuver. Optimal regulator theory (outlined 
above) was applied in an attempt to reduce the extreme acceleration peak. In princi- 
ple, the theory worked very well. However, the reduction in lateral acceleration 
during the first 3 to 5 seconds of the rolling maneuver was almost inconsequential. 
Extreme variations in the weighting matrices, Q and R t , of the performance index 
railed to change the results. 


An open- loop control input was applied to alleviate the lateral acceleration. 
The aileron was deflected, at maximum rate, to its stop (30°), held for awhile, then 
returned to near zero deflection angle. From consideration of equation (2), an ini- 
tial positive rudder deflection is required. The response to such control inputs is 
shown in figure 6 . The control inputs are shown in figure 6 (b). The airplane rolled 
u 4 , , the re< l uired 3 seconds (See Ref. 4) and the lateral acceleration peak at 
t e p ot s station was less than .075 g. Also, the peak occurred at about 2 seconds 
rather than immediately after the control input. However, lateral acceleration 
changed signs after 3 seconds. The responses of yawing velocity and sideslip angle 
however, were not satisfactory. Sideslip angle increased to about 15° and yawing 
velocity was negative for 3 seconds. 

The undesirable (ay)pg seems to occur in the command response because of the 
sudden control input. Therefore, crossfeed and prefiltering, to slow up the input 
and coordinate the controls, seems to be required. Figure 6 illustrates the desir- 
able crossfeed effects • 


Two aileron and rudder interconnects were placed in the forward control path. 
The gain matrices were chosen to improve roll and yaw control by negating rolling 
moment caused by rudder deflections and by negating yawing moment caused by aileron 
deflections. The feedforward matrices are 


and 




-.2538 

1.0 


Aileron deflections produced no yawing moment when Gi was used, 
rudder deflections produced no rolling moment when G 2 was used. 


In addition. 


Response of the airplane to the control input of figure 6 (b) with Gi and Go 
in the forward control path is shown in figure 7. Only minor improvements in the 
response were obtained by using G]^ or G 2 (compare Fig. 6 with Fig. 7). Further- 
m °re, neither G^ nor G 2 was superior. 


CONCLUDING REMA RKS 


The resuits of this study show that modal control theory provides desired stabi- 
lity and desired transient response characteristics. Proposed criteria for the 
maximum allowable lateral acceleration at the pilot station was satisfied by using 
open-loop control inputs. However, large adverse sideslip angle excursions were 
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required; and the yawing velocity was negative for about 2 seconds. Optimal regula- 
tor theory did not give a significant reduction of the lateral acceleration level 
during transient motion. 


APPENDIX 


Theory 


The problem of selection of characteristic values 
briefly outlined as follows (See Refs. 7 and 8). 


and characteristic vectors 


is 


Consider the controllable system given by equation (1). That is. 


x = Ax + Bu 

where xeR^, ueR 2 , and A and B are properly dimensioned. Assume that B is full 
rank. The state variable feedback law takes the form 


u = Kx + u 
c 


(A-l) 


where K is the feedback matrix, 
system matrix A+BK satisfies 


The problem is to select K so that the closed- loop 


[A+BK] v ± - X^ (i=l ,2,3,4) 


(A-2 ) 


where ^is the ith characteristic value and Vi is the corresponding character- 
1st lc vector • 


Partition equation (A-2) as 


'n 


12 


A • A 
21 A 22 


B„ 


fv K 2 ] 


= A. 


w. 


( i = l >2,3,4) 


(A-3) 


where An, Bi, and % are 2x2 matrices and other matrices are compatibly dimen- 
sioned. Also, note that Bi, is nonsingular. z t is a 2x1 vector, z, and w, 

are partitions of the characteristic vector, v i , and v< T « rz.« T *v..Ti 

Completing the multiplication of equation (A-3) yields 1 1 1 ’ 1 


[ A U + hh K + t A 12 + hhK " X i z i 

[A 2 i + B 2 K l ]« 1 + [A 22 + B 2 K 2 ] Wi - X lWl 


(A-4) 

(A-5) 


Solve equation (A-4) for Ki Zi + K 2 wi and substitute into 
obtain the following constraining relationship 


equation (A-5) to 


fVn-m ‘ F ^ w i " f G + V^ z i U-1,2,3,4) 


(A-6) 


12 


where I n _ m *l 2 is a 2nd order identity matrix and S, G, 
by 


and F are matrices defined 


S ■ V? 


G - A 21 - SA 
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(A-7) 


F 


A 22 ~ SA 12 


Equation (A- 6 ) constitutes a set of n-m=2 linear equations in the n=4 unknown 
elements of each characteristic vector. Thus, if \ ± is not a characteristic value 
of F, at most m elements corresponding to the z ± vector can be arbitrarily chosen. 
Then the remaining n-m characteristic vector elements corresponding to the w< 
vector can be computed from equation (A- 6 ) as 


C i z i 


(i = l *2,3,4) 


(A- 8 ) 


where 0 ^= [ X£l n _ m -F ] ^[g+X^s] is defined as the "modal coupling matrix" 
corresponding to X ± . f b 

The preceding analysis shows that all n characteristic values and up to n*m 
characteristic vector elements can be arbitrarily assigned through the use of state 
variable feedback. The resulting modal matrix 


[W-" :v n] 


however, must be nonsingular. 

With a nonsingular modal matrix, V, and characteristic value matrix. A, chosen 
subject to the constraints of equation (A- 6 ), the closed-loop system matrix A=*A+BK is 
uniquely determined by 


*11 


A 

: ^2 


21 


The matrix A is given by 


22 


VAV 


-1 


(A-9) 


a 

-P 

0 


a 


4 -I 


(A-10) 


where Xi and X 4 are real and X 2 =ort-jB and X 3 »a-j 8 . 

The required feedback matrix, K, which yields this closed-loop matrix, A, can be 
easily computed by using the relations 
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( A— 1 1 ) 


*1 ' B i' l [*n - V 

k 2 - Bl - l [A l2 - a 12 ] 


K - [k x :K 2 ] 


(A-12) 
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TABLE I.- MASS AND DIMENSIONAL CHARACTERISTICS OF 
SUPERSONIC TRANSPORT AIRPLANE 


Weight, N 1 924 479 

Reference wing area, m 7 g 4 75 

Wing span m 38 : 6 6 

Wing leading-edge sweep, deg (see Fig. 1) 74.00/70.84/60.00 

Reference mean aerodynamic chord, m 27.00 

Center— of— gravity location, percent c ......... 55 

Static margin, percent * - 3.9 

I X» kg " m * * * 6 887 550 

I Y’ kg-m 67 994 260 

2 

1 V kg-m 72 902 230 

2 

I XZ’ kg-m -2 833 660 

Maximum control surface deflections: 

6 f» deg 0 to 40 

6 a» deg ±30 

6 af» deg ±22.5 

V deg * • ‘ ±50 

6 r’ deg ±35 

Maximum control surface deflection rates: 

S f , deg/ sec 


6 , deg/sec . 

3 • • • • 

, deg/ sec 

6 , deg/sec 

g ******* 

6 r , deg/ sec 

Horizontal tail: 

Gross horizontal-tail area, 

Mean aerodynamic chord, 

Distance from center-of-gravity to horizontal-tail 0.25c, m 
Vertical tail: 

Exposed vertical-tail area, m 2 

Mean aerodynamic chord, m 

Distance from center-of-gravity to vertical-tail 0.25c, m . 


. ±70 
. ±40 
. ±50 
. ±50 

49.80 

6.04 

32.90 

16.72 

6.35 

36.41 
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TABLE II.- AERODYNAMIC CHARACTERISTICS OF 
SUPERSONIC TRANSPORT AIRPLANE 


Airspeed - 78.71 m/sec; Altitude = 91.44 m; a - 8°; 6 f - 40° 


o 

II 

1.1793 rad*" 1 

C Y " 
r 

.4154 rad" 1 

c - -.00723 deg" 1 

8 

C A - 
p 

-.1389 rad" 1 

c * - 
r 

.1946 rad" 1 

C. = -.00219 deg" 1 

C 

n 

P 

-.0747 rad" 1 

C 

n 

r 

-.2941 rad -1 

C « .00160 deg -1 

P 

II 

u 

c o 

>< 

o 

.00100 deg -1 




c l = 
*6r 

.00017 deg -1 




C 

n 6r 

-.00119 deg -1 




C Y = 

6s 

—.00008 deg 1 

C Y = 

6af 

—.00026 deg 1 

C = -.00026 deg -1 

X «a 

C* " 
6 S 

.00012 deg 1 

C* 

6af 

.00040 deg 1 

C, = .00036 deg 1 

6a 

C 

n 6s 

0 

C 

n 6af 

.00009 deg -1 

C = .00014 deg 1 

n 6a 

NOTE: 

The aileron, spoiler, and 

r.i i 

flaperon effectiveness terms were combined 


follows : 


'() 


1.667 C 


6a 


() 


+ .75 C 


6s 


() 


+ C 


6af 


() 


6a 


since 6 g -1.6676 a and 6 a j=0.756 a for this report. 
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TABLE III.- COMPARISON OF AIRPLANE STABILITY AND REQUIRED STABILITY 


Roll Mode: 


Spiral Mode: 


Dutch-Roll Mode: 


♦Whichever gives 


Airplane Stability Required Stability 


5* 

II 

1.64 

t r < 

1.4 

T l/2 = 

22.4 seconds 

T 2 > 

20 seconds 

C d = 

.093 

? d > 

.08* 

W nd = 

.825 

%d > 

.4 

d“nd = 

.077 

5 d“nd > 

.15* 


the largest value of £ . 

d 


17 




46.37 


Wing reference plane 

















40 


30 


0 . 

deg 


20 


10 


N3 

Ui 


0 


10 


P. 

deg/sec 


5 


0 


2 


e, sec 


(a) State response. 


Time, sec 


5.- Response of augmented airplane to an aileron input. 




Time, sec 


Time, sec 


(b) Lateral acceleration and control deflections. 


Figure 5.- Concluded. 








40 


Feed forward matrix 



Time, sec 


Time, sec 

(b) Lateral acceleration and control deflections. 
Figure 7.- Concluded. 
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